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Abstract 

The emission of vv pairs off electrons in a polarized ultra-intense electromagnetic (e.g. laser) 
wave field is analyzed. We elaborate on the significance of non-linear electrodynamics effects (i.e., 
multi-photon processes) and the peculiarities of neutrino production. Special attention is devoted 
to the convergence of the reaction probabilities as a function of the number of absorbed photons. 
Expressions for large field intensities are provided. The asymmetry between the probabilities of 
electron and {jl+t neutrino production depends on initial conditions such as energy of the wave field 
photons and the field intensity. These findings differ from the lowest order perturbative calculation 
of the reaction 7 + e — > e' + vv. 

PACS numbers: 13.35.Bv, 13.40.Ks, 14.60.Ef 
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I. INTRODUCTION 



The exact solution of Dirac's equation for an electron moving in the field of a plane 

n 

electromagnetic wave was found by D. M. Volkov in 1935 The electron wave function, 
compared to the field-free case, changes due to a modification of its spinor structure and 
the appearance of an additional phase factor. The electron momentum changes to an effec- 
tive quasi- momentum, and the electron mass becomes an effective "dressed" mass. These 
modifications depend on the dimensionless variable £ 2 related to the amplitude of the elec- 
tromagnetic four-potential 

e = -^!> (i) 



e 

2 



where e is the absolute value of electron charge (e = Ana with a 



the electron mass. (For a manifestly gauge invariant formulation cf. {3].) 

Decades later, Volkov's solution was applied to Compton scattering 
positron pair production [3] in strong electromagnetic fields. A consistent systematic anal- 
ysis of these electromagnetic and further weak processes, such as pion and muon decays, 
vv emission by an electron in an external field etc. was performed by Nikishov and Ritus 
and coworkers in a series of papers js il] and summarized in the review l^ . 



/137.035) and M e is 



-|6| and electron- 



Later, some 



aspects of weak interaction, in particular neutrino emission, by electrons in a strong elec- 



tromagnetic field were considered in Refs. 13-151. The twofold extension of QED for strong 

n 

electromagnetic fields was discussed in the recent paper [16j. 

The main result of these previous studies is the conclusion that the quantum processes are 
modified significantly in strong electromagnetic fields. For instance, an electron can absorb 
or emit simultaneously a certain number of field photons, depending on the initial conditions 
of the considered process. This fact, together with the modifications of above mentioned 
electron properties, results in strong non-linear and non-perturbative effects which can not 
be described within the usual perturbative quantum electrodynamics (pQED). Consider, 
for example, the emission of a photon with four-momentum k' by an electron moving in 
a electromagnetic wave field. The process depends on the invariant variable u = |^ 
which varies in the range of < u < u n = j^"^!^) f° r the absorption of n photons with 
four-momenta k ~ (ul, kt,) by the electron with four-momenta p ~ (E e , p e ) and p 1 prior 
and after the emission process. One can see that (i) the kinematical limit u n (phase space) 
increases with the number of absorbed photons ("cumulative effect") and (ii) decreases with 



increasing field intensity £ 2 because of the electron mass modification. On the other hand, 
the contribution of higher harmonics also increases with £ 2 , where, following j^j], we use the 
notion "harmonics" for processes with different n's. 

Since £ 2 plays an important role, it seems to be useful to recall the relation between £ 2 
and the electromagnetic (laser) field intensity /, where the electromagnetic field is considered 
as a classical background field. For the case of a monochromatic circularly polarized plane 
wave with four-potential = (0, A), where A(0) = a x cos<p + a^sin^, <fi = k ■ x, and 



&5 a^a^ 



0, i.e. k ■ x = uit — kxx, the average value of A 2 is equal to 



meaning £ 2 



On the other hand, the field intensity may be expressed through the 



electric (E) and magnetic (H) field strengths by / 



fE 2 + H 2 



cE 2 . Taking into account 



E = —dA/dt, one gets an expression for the average intensity / in terms of the amplitude 
a, I = ca 2 u 2 , which leads to £ 2 = t ^°^ 2( , A| /, where \l = 2nhc/ujL is the wave length of 
the electromagnetic field A^ 1 . The dependence of £ 2 on the electromagnetic field intensity / 
for different wavelengths Al is exhibited in Fig. [TJ 
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FIG. 1: Dependence of £ 2 on the electromagnetic wave field intensity / for three values of wave 
length Xl- 



The wavelength = 0.8 fim (or 



1.55 eV) corresponds to the widely used titanium- 
sapphire laser oscillator (cf. Refs. 17H20|). The short wavelength A^ = 2.5 nm (or ojl — 
0.5 keV) corresponds to the soft x-ray (SXR) free electron laser at SLAC [211] . The long 
wavelength Al = 40 /im (or oj j ^ 0.03 eV) may be obtained at the free electron laser for 
infrared experiments (FELIX) 22[. One can see that £ 2 varies within a fairly large range, 
depending on the field intensity and wavelength. 

In the low-frequency limit, ujl — > 0, the intensity parameter becomes large, i.e. £ 2 — > 
oo at fixed intensity / or E. This limit was considered in some detail by Nikishov and 
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Ritus 



BE 



12] who pointed out that the invariant variable 



remains finite and the total probabilities of most of the considered processes depend only 
on x Ql ■ Here, F^ v = d u A^ — d^A v is the electromagnetic field tensor. Such a case of 
simultaneous limits of £ — > oo and w^, — )■ at finite I corresponds to the situation of an 
electron interacting with a constant (crossed) electromagnetic field. 

Note that two asymptotic regions of the external field were considered in most of the 
above quoted papers. One corresponds to the weak-field limit ^ 2 < 1. In this case, only 
a limited number of harmonics n < 2 contributes. The opposite case of large intensity 
£ 2 — > oo with oul — > allows for two asymptotic limits: x ^ 1 an d X ^ 1- Of course, such 
an analysis of limiting cases is interesting and important by its own. However, the rapidly 



evolving laser technology 23| can provide conditions where the limit of £ 2 ^> 1 is achieved 
at finite ul, as well as x ~ 1 as can be inferred from Fig. [Hand by numerical evaluation 
of Eq. (j2J). Therefore, it seems relevant to consider the probabilities of quantum processes 



without the restrictions imposed in 12H15] 



The goal of present work is accordingly an analysis of neutrino pair emission off an 
electron moving in a strong external electromagnetic (laser) wave field in a wide region of 
£ and x- O ur paper is organized as follows. In Sect. II, we consider the neutrino pair 
emission. A scheme is presented to overcome convergence problems in the expansion in 
terms of harmonics. The employed method is similar to the one-photon emission process 
which is outlined in Appendix A to expose these similarities and the important differences. 
The perturbative neutrino pair emission is recapitulated in Appendix B. Our conclusions 
can be found in Sect. III. 

II. EMISSION OF A NEUTRINO PAIR 
A. Basic formulas 

Similar to the emission of a photon by an electron moving in an electromagnetic external 
(background) field (see Appendix IA II) one can evaluate the emission of a neutrino pair ViVi 
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e(q) e'(q') 



FIG. 2: Diagram for the emission of a neutrino pair off an electron in an external wave field with 
effective low-energy eevv vertex. The double lines depict electron Volkov states. 



of species i = e, fi,r by the S matrix element 



qii) 



V2 



d A x 



(3) 



y / 2E u 2Ep 

corresponding to the diagram in Fig. [2j Mq is the invariant mass of the vv pair: Mq = 
Q 2 — {.Vu + Vv) 2 with p u ~ (E u ,p u ) and p„ ~ (E p ,p p ) as four-momenta of neutrino and 
antineutrino. (We employ units with h = c = 1.) The Volkov solution ipi(x) in ([3]) (cf. 
Eq. f]A2|) in Appendix IA ip describes the in-state of the electron in an external wave field, 
while ip*f{x) refers to the out-state of the electron also accounting for the external field. The 
neutral neutrino current 



(4) 



couples directly to the electron current ^j7 a (Cy' ) — C^jsjipi with a strength given by 
Fermi's constant Gf = 1-66 • 1CT 5 GeV~ 2 . The expression ([3]) holds in the local limit where 
all momenta involved in the process are much smaller than the masses of the intermediate 
vector bosons Z° and W . Then, one can obtain "universal" interactions described by the 



effective low-energy Lagrangian 



24| for the direct current- current interaction 



C 



(i) 
eff 



Gf 

71 



(5) 



with 



C 



c 



(<') 



- + 2 sin 2 6 W , 
1 



C 



c 



v 



- + 2 sin 2 9 W 
1 

"2 ; 



(6) 
(7) 



the average value of sin 2 9w ~ 0.23 is taken from Ref. 25|- In the weak-field approximation, 
where the interaction with the external wave field is mediated by one photon, the diagram 
in Fig. [2] would be resolved by diagrams exhibited in Fig. [3] corresponding to the process 



7 + e — > e' + vv which can be dealt with perturbatively (see Appendix [B]) . The effective weak 
vertices for the eevv interaction, in turn, are resolved in the tree-level approximation within 
the standard model as exhibited in Fig. HJ The advantage of the diagram in Fig. [2] and 
Eq. ([3]) is that multi-photon effects are included which become important at high intensities 
where the external wave behaves more and more as a classical field. 




(a) (b) 

FIG. 3: Lowest order diagrams for the reaction 7 + e — > e' + vv with effective low-energy eevv 
vertices. Diagrams (a) and (b) correspond to the direct charge and neutral vector boson exchanges. 
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FIG. 4: Lowest order diagram for the reaction 7 + e — > e' + vv. Diagrams (a) and (b) correspond 
to the direct charge (neutral) vector boson exchange; (c) and (d) describe corresponding crossed 
channels. 

As in the case of the emission of one photon (see Appendix [A]), the differential probability 

of neutrino pair emission can be represented by the infinite sum of the partial contributions, 

00 

dW = J2 dW {n) (8) 

n=l 

with 

2(8tt) 3 q (l + u) 2 2vr ' V ' 

R (n) = l(Q a Qp-g a pQ 2 )Mf , (10) 
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where <p q i is the azimuthal angle of the outgoing electron and 

Analog to the case of one-photon emission (see Appendix IA ip the label n refers to the 
number of photons absorbed from the external field. For the sake of simplicity, we skip 
henceforth the index %. The electron tensor 

Mf = Tr[(/ + M e )S a n (C v - C Al5 )(j> + M e )S^C v - C Al5 )] (12) 

incorporates the electromagnetic wave field via 

bn ~V 2k.pk.p') Bn + ^ Me [\2k-p' + 2k.p) Bn + \2k.p' + 2k-p) Bn _ ■ 
Here, n^) is the unit vector of ai(2), and the functions B n are related to the Bessel functions 
(cf. Refs. flQ): 

= J n (z)e in ^' , 

B< 2 > = I (j fl+1 e i <»+ 1 >«V - J^e* 1-1 **) . (13) 
The argument of these Bessel functions is 



' « ( ! - « )- I± ^, , M L, , (14) 



where u n is the kinematical limit of the invariant variable u defined by 

_ 2n(k-p) 

n "M e 2(l + a { ] 

which determines the upper limit of Mq at given u by 

mi =M e 2 (i+e 2 )^^. (i6) 



The functions By^ in Eq. f JTOj) do not depend on </> 9 /; the explicit expression reads 

= C 2 V F V + C\F A + 2\C v C A Fj, (17) 

F v = -M 2 (2M 2 e + M 2 Q )J 2 n + e 2 MjMl (l + ^^y) AJ 2 , (18) 

F A = Mq (AM 2 - M 2 )J 2 n + i 2 M 2 E ( Mq + (Mq + 2M 2 )—^-—] AJ 2 n , (19) 



2(l + u 



M'A I 1 - ; M " ; ] - M, 2 



, ^(2 + ^) V ^ 2(u n -u 

7(1 +«) 



Ff = M 2 M e ^= ^ ^ ^ -Lf- '-J n (J n+1 - J n _ x ) , (20) 



V max V , 



where A = ±1 is the relative phase (polarization) of the amplitudes 
Jn+i + Jn-i~ 2 J%. Our expression coincides with the result of Ref. 



ai and &2, and A J 2 
15| in a different notation. 



B. Numerical evaluation 
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FIG. 5: Differential probability of neutrino pair emission as a function of logu for £ 2 = 0.1, 1 
and 10, shown in left, middle and right panels, respectively, for the lowest harmonics n < 4. The 
result for the reaction 7 + e — > e + vv calculated perturbatively (see Appendix B) is shown by the 
dashed curves labeled by "free". The numbers (from 1 to 4) correspond to the number of absorbed 
photons. The energies of the external (laser) photons and incoming electrons are chosen to be 
ujl = 1.55 eV and E e = 40 MeV, respectively. The invariant mass of the outgoing neutrino pair is 
fixed at Mq = 10 eV. 



In Fig. [5] we exhibit the differential emission probabilities summed over all neutrino types 
for a head-on collision of 40 MeV electrons with a laser beam characterized by Xl = 0.8 fim 
(or ujl ~ 1.55 MeV). We fix the invariant mass of the outgoing neutrino pair by Mq = 10 eV. 
The numerical results for the first two harmonics coincide for £ 2 = 0.1 with the prediction 
15| based on the asymptotic decomposition of W(u) for n — 1, 2 in the limit £ 2 < 1.) The 
prediction of perturbative QED for the reaction 7 + e — > e' + vv (see Fig. |3] and Eq. (IB5j) in 
Appendix |B]) is shown by the dashed curves. One can see that for small £ 2 the result for n = 1 
coincides practically with that of pQED. Similar to the case of non-linear Compton scattering 
(see Appendix I A ip . even in this case the contribution of higher harmonics increases the phase 
space, i.e. the kinematic limit, and modifies the total probability. When £ 2 increases, the 
difference between the non-perturbative calculations and pQED is rather large, even for 
n — 1, see middle and right panels of Fig. |5j 
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FIG. 6: The total probability of neutrino pair emission as a function of £ 2 for the first 25 harmonics. 
The thick curve is their total contribution. The electron and photon energies are chosen as E e = 
40 MeV and ujl = 1-55 eV, respectively. 

In Fig. E]we show the total probability normalized as Wj p e , where p e is the initial electron 
density divided by the electron mass in the electron rest frame, 



for the sum of all types of neutrinos as a function of £ 2 for the first 25 harmonics. For 
small £ 2 , £ 2 < 0.1, along with the predominant contribution of the first two harmonics the 
contribution of the higher harmonics is significant. When £ 2 increases, e.g. for £ 2 ~ 10, 
the contribution of higher harmonics exceeds the contribution of lowest harmonic by orders 
of magnitude. The qualitative difference in the relative contribution of higher harmonics 
to the total probability of the emission of a photon (cf. Fig. [13] in Appendix A 1) and 
neutrino pairs is explained by the employed four-fermion structure of the weak-interaction 
eevv vertex and the three-particle phase space. Therefore, the problem of convergence 
for the total probability in case of neutrino emission becomes severe and deserves special 
consideration. 




(21) 
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C. Overcoming convergence problems 



Figure [7] illustrates the convergence problem of the total vv emission probability with 
increasing £ 2 for the case of ul = 1.55 eV and E e = 40 MeV, where the number of included 
harmonics goes up to n max = 140. One can see some saturation at £ 2 < 2. However, for 
£ 2 > 10 the difference of probabilities with n max = 25 and n max = 140 is more than two 
orders of magnitude. 
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FIG. 7: The total probability of neutrino pair emission as a function of £ 2 for different values of 
n max from 25 to 140. Kinematics as in Fig. [6l 



The problem of convergence may be solved by the same method as proved useful for 
the one-photon emission (see Appendix IA 11) using the transformations in Eqs. (1A10|) and 
(lAllj) . The main difference to the one-photon emission (cf. Eq. ( 1A13I) ) is a modification of 
the argument of the Airy functions, which emerge from the above Bessel functions, as 



y = t 1 + r z + 



1 + wM 2 



(22) 



u 2 M 2 , 

and an additional integration over the invariant mass g?Mq. Also, the interference term in 
Eq. ( 120]) proportional to the rapidly oscillating combination $ (?/)$'(?/) is negligible and can 
be omitted. The final expression for the probability of vv emission in the limit of 

n mai — > oo reads 



48tt 5 



y/idu 



dM, 



Q 



dr ( Fy^Cy + F^Cl 



(23) 



10 



4 A) = -Mg(2M 2 + M 2 )d> 2 (y) + \m%MI (\ + ^^y) (y^ 2 (y) + ^(y)) , 

= Mj(4M 2 - M 2 )$ 2 (y) + ^M 2 ^Mg + (Mg + 2M 2 )^-^ (y$ 2 (y) + $' 2 (y)) , 
where t = (u/2x) 2 ^ 3 - 
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FIG. 8: The total probability W^ A \i, \) of neutrino pair emission as a function of \ f° r different 
values of £ 2 as indicated in the legend. The stars and crosses correspond to the asymptotic values 
of W(oo, x) calculated by Eq. (f2i|) for % ^ 1 an d ^ respectively. 

The emission probability W^ A \^,x) for all three types of neutrinos summed up is shown 
in Fig. |H] for a wide region of x an d £• The stars and crosses correspond to the asymptotic 
values of W(oo, x) for \' < 1 and x>l, respectively |l4j : 

' P -^f(f(C 2 v + C\) + 63Cl), X «l, 

(24) 



W(oo,x) = < 



216 



^L(C 2 +Ci) (In x - 0.577 -§ In 3 -|), x»l ■ 



Now, the dependence W^ A \^,x) on £ is weaker compared with the photon emission (cf. 
Fig. [To]) , mainly because of the additional integration over Mq. Note that the asymptotic 
estimates by Eq. ( |24p do not match in the intermediate region 0.1 < x ^ 30- 

In Fig. [9] we show the total probability W(£,x) of neutrino pair emission for all neutri- 
nos calculated for wide initial experimental conditions ranging from ujl = 0.03 eV up to 
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FIG. 9: The summed total probability of emission of all three types of neutrinos as a function of £ 2 
for a finite number of harmonics, n max = 25 and 140 (dashed curves), together with the asymptotic 
probability given by Eq. (|23p (solid curves). The electron energy is chosen as E e = 40 MeV, and 
the laser photon energies ul are 0.03 eV (left panel), 1.55 eV (middle panel) and 0.5 keV (right 
panel). 

0.5 keV as a function of £ 2 , where £ 2 varies from 0.1 up to 10 3 . This interval covers possible 
experimental conditions illustrated in Fig. [TJ At £ 2 > 6, W(£, x) is evaluated using the 
asymptotic expression of Eq. fl23|) . For £ 2 < 6, the probability might be evaluated as a sum 
of partial harmonics, in our case up to n max = 140. One can see that, at large £ 2 , the differ- 
ence between the probability calculated as a sum of a large but finite number of harmonics 
and its asymptotic value which includes an infinite number of harmonics is several orders of 
magnitude. 

Finally, we would like to note the following. In average, the difference of neutrino pair 
emission probabilities for ojl = 0.03 eV and 0.5 keV is about 17 orders of magnitude (cf. 
left and right panels in Fig. [9]). This is much greater than the corresponding difference in 
one-photon emission shown in the left and right panels of Fig. [16], where the corresponding 
difference is about 5 orders of magnitude. The difference between the two processes is 
explained by the different x dependence of W(£,x), shown in Figs. Eland [HI respectively. 
The average values of x for ul = 0.03 eV and 0.5 keV are 10 -4 and 4, respectively. Therefore, 
W(C,, x) increases with u^. The sharp increase of the probability with x i n cas e of neutrino 
pair emission is a consequence of the strong energy dependence of the total probability which 
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can be traced back to the four-fermion structure of the eevv matrix element. 



D. Asymmetry of v e v e vs. Vut^ut emission 
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FIG. 10: Asymmetry Ai e ^ T \ as a function of £ 2 for different values of the laser photon energy u>l- 
Left panel: low £ 2 ; probabilities are calculated as a sum of a finite number of harmonics according 
to Eq. 0. The prediction for the 7 + e — > e' + vv reaction is shown by the dot-dashed curve 
labelled by "free". Right panel: high £ 2 ; the probabilities take into account an infinite number of 
partial harmonics according to Eq. 



It seems to be interesting to analyze the asymmetry of the emission of electron and muon 
plus tau neutrino pairs which may be defined as 

W^^W). (25) 

Corresponding predictions for two regions of £ 2 are shown in Fig. [10] for different values of 
the wave field photon energies. In the region of £ 2 < 1.2, the probabilities are calculated 
as a finite sum of partial harmonics with n max = 50 (cf. Eq. (J2J)). The prediction for the 
7 + e — > e' + vv reaction (see Appendix |B]) is shown by the dot-dashed curve labelled by 
"free". Predictions for ui = 0.03 and 1.55 eV are virtually identical. For £ 2 1, all 
predictions are close to the results for the 7 + e — > e' + vv reaction: A^t) — —0.07. All 
asymmetries are small in absolute value and negative, meaning a slight dominance of the 
emission of \i + r neutrinos compared to electron neutrinos. 

At large £ 2 , the probabilities are calculated as an infinite number of partial harmonics 
(cf. Eq. fl23D). For ul = 0.03 and 1.55 eV the range of variation of \ is 2-10~ 5 < \ < 3-3-10" 4 
and 1.6 ■ 10~ 3 < x ^ 3 • 1.5 ■ 10~ 2 , respectively. That is, the variable \ is small and 
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therefore, the asymmetries at these two energies are close to each other and close to the value 
-4( e ,/ir) — —0.17 which results from the asymptotic expression Eq. f l24l) . For lol = 0.5 keV 
and 0.34 < \ ^ 4.8 the asymmetry is positive and increases with £ 2 . Note that, within the 
considered range of £ 2 , Eq. (I24"|) does not apply and A( e ^ T ) is smaller than its asymptotic 
value ~ 0.4 which would be found from Eq. f)24p . We predict a distinct dominance of 
emission of electron neutrino pairs compared to the sum of muon and tau neutrino pairs in 
the keV-range of the photon energy and at large £ 2 . The asymmetry increases as ln£ 2 in 
the interval of 10 < £ 2 < 1000. 

Finally, we note that our non-perturbative calculation of A( e ,nT) is strongly different from 
the prediction of pQED (cf. dot-dashed curve in the left panel of Fig. [TD] ) in all considered 
intervals of ojl and £ 2 unless £ 2 1. 

III. SUMMARY 

In summary we have considered z/jZ/j (i = e,fi,r) emission off an electron in a strong 
electromagnetic wave field in a wide range of energy of the wave field photon energy ul and 
reduced field intensity £ 2 . Similarly to previous work we expressed the emission probability 
as a sum of partial harmonics, where each harmonic describes the interaction of an electron 
in-state with n field photons coherently. For the first time, we made a summation over 
harmonics up to a quite large number of harmonics and found that, at large values of 
field intensity, £ 2 > 10, which can be achieved in current and future laser facilities, the 
convergence is rather weak. Therefore, we have elaborated a method allowing for a complete 
summation of all partial harmonics. The method is tested for one-photon emission, i.e. non- 
linear Compton scattering. Using this new approach we calculated neutrino pair emission in 
a region of ul and £ 2 which can be reached experimentally in near future. We have shown 
that, at large £ 2 , the difference between the finite and complete sums of partial harmonics 
reaches a few orders of magnitudes. 

In case of neutrino pair emission we also analyzed the non-trivial asymmetry between 
the production of electron and fi + r neutrino pairs. We found that the asymmetry depends 
strongly on initial conditions expressed via ul and £ 2 . At low ul, the asymmetry is negative 
corresponding to the dominance of emission of the fi + r neutrinos, while at large ojl the 
asymmetry changes the sign indicating the dominance of electron neutrino pairs. 
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Finally, we note that all calculations (and conclusions) have been done for the sake of 
simplicity for the initial electron energy E e = 40 MeV, which corresponds to the energy 
of the superconducting electron accelerator ELBE in FZ Dresden- Rossendorf 26|. One of 
the key variable in the considered processes, x = £ k ' P/M^, where k and p are the photon 
and electron four- momenta, respectively, directly depends on E e . Therefore, the values 
of emission probabilities would also depend on E e . The corresponding evaluation of this 
dependence as well as the analysis of other processes will be done in forthcoming papers. A 
further step towards realistic estimates is related to the inclusion of a temporal shape of the 



external laser wave field, as considered e.g. in 



23-|29|. 
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Appendix A: Photon emission off an electron in a strong electromagnetic wave field 



The methods employed in Section II are guided by the one-photon emission process off a 
Volkov electron, i.e. the non-linear Compton effect. To expose the similarities and differences 
we recall the essential steps and clarify further the notation. 

1. Strong external field 

Y'(k')^< 
e(q) e'(q') 

(a) 

FIG. 11: (a) Diagrammatic representation of the emission of a photon off an electron in an external 
wave field. The double lines depict Volkov states, (b) and (c) are direct and exchange terms of the 
lowest-order diagrams for the perturbative treatment of the Compton process. 

Let us consider the emission of a photon off an electron moving in a plane electromagnetic 
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wave described in the Introduction. This process is described by the S matrix element 

Sfi = -ie / ^(7 • 4)^*4= , ( A 



Jk'-x 

where u' is the energy of the emitted photon with four-momentum k' and ipim is the electron 
wave function in the initial (final) state given by Volkov's solution of the Dirac equation 

k-x 



1 + 



2(k-p) 



exp 



e(p ■ A) 
WW 



u 



V Q -iq-x 



where 7 denote Dirac matrices and q is the quasi-momentum 

e 2 (A 2 ) 



e 2 a 2 



2{k-p) * 2{k-p) 
of the dressed electron with effective mass 

e 2 (A 2 ) 
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Note that Eq. (IA1I) employs the Furry picture: The field A^ is considered as external classical 
(background) field, and the emission of a photon with wave four-momentum k' and polar- 
ization Ef is described in lowest order of perturbation theory, see diagram (a) in Fig. [TTJ 

The dependence of the potential A on k ■ x in Eq. (1A2|) results in the following structure 
of the S matrix element 



S fi 



M fi (kx)e 



h v/2g 2^2u/ 
where a Fourier decomposition is used: 



i(q-q'-k>) 



■ x d 4 x= M 1 } i {2Ti) A 5 A {q + nk - q' - k'),(A5) 



M(kx,k,k',q,q') = e ~ ink ' X M n (k,k',q,q') 



(A6) 



Thus, the amplitude is represented as a sum of an infinite number of terms which are referred 
to as partial harmonics. Each harmonic can be attributed to the absorption (emission) of 
n photons from (into) the external field A characterized by the wave four-vector k. For the 
photon emission off an electron the corresponding conservation law reads q + nk = q' + k', cf. 
Eq. (1A5|) . In case of on-shell photons, k 2 = k' 2 = and k' ■ q' > hold and therefore n > 1. 
Correspondingly, the differential probability is the infinite sum of partial contributions: 



dw = J2 dw(n) > 



(A7) 



n=l 



1(3 



where the partial harmonics W^ 1 ' are expressed through Bessel functions J n of the first 
kind Q, 3| 

dwin) = ^ (TT^F {- 2J * {z) + + 2(TT^) ) ( J "+ lW + J «- l( ^ " 2J « ( ^ } (A8) 



with 
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k ■ p 1 E' e + u' cos 9 
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FIG. 12: Differential probability of one-photon emission as a function of log it for £ 2 = 0.1, 1 and 
10, shown in left, middle and right panels, respectively, for lowest harmonics n < 4. The pQCD 
Compton scattering is shown by the dashed curves labelled by "K-N". The numbers (from 1 to 4) 
correspond to the number of absorbed photons. The energies of the laser photons and incoming 
electrons are chosen to be ojl = 1-55 eV and E e = 40 MeV, respectively. 

As an example, we exhibit in Fig. [12] the differential probabilities for photon emission by 
40 MeV electrons colliding head-on with a = 0.8 /im (ujl — 1.55 MeV) laser beam. This 
process can be studied experimentally, e.g., at superconducting electron accelerator ELBE 
in conjunction with the 150 TW laser Draco in FZ Dresden- Rossendorf 26]. The solid 
curves correspond to the partial contributions which come from the coherent absorption of n 
photons with n = 1...4. The prediction of perturbative QED Compton scattering, described 
by the Klein-Nishina formula (see Appendix IA2I) . is shown by the dashed curves. One can 
see that for small £ 2 , the results for n = 1 coincide practically with the result of pQED, 
cf. also However, even in this case the contribution of higher harmonics increases the 
phase space and modifies the total probability. When £ 2 increases, the difference between 
non-perturbative calculations and pQED is rather large, even for n — 1. In general, the 
modification of the kinematical limit follows Eq. (fl~5|) : u n increases with n and decreases 
with £ 2 . 
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In Fig. [13] we show the total probability of one-photon emission as a function of £ 2 for 
the first 25 harmonics. For convenience, we show reduced probability W/p e , where p e is the 
density of initial electrons divided by the electron mass in the electron rest frame. 



> 
CD 
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FIG. 13: The normalized probability of one-photon emission as a function of £ 2 for the first 25 
harmonics. The thick curve depicts the completely summed up result. The electron and photon 
energies are chosen as E e = 40 MeV and ujl = 1-55 eV, respectively. 



On can see that for small £ 2 , £ 2 <C 1, the total probability is saturated by the contribu- 
tion of the first two harmonics. However, at large £ 2 , £ 2 > 10, the contribution of higher 
harmonics becomes large and the convergence of the total probability as a function of the 
number of harmonics n is weak. Figure [TH illustrates the convergence of the emission prob- 
ability with increasing £ 2 for the case of ujl = 1.55 eV and E e = 40 MeV, where the number 
of harmonics increases up to n max = 140. One can see some saturation as long as £ 2 < 10. 
When £ 2 increases, n max must further increase, in principle, as n max ~ £ 3 [l^ . 

For the calculation of the total probability for large but finite £ 2 we use the method of 
Ref. 12j based on utilizing the properties of the Bessel functions J n (z) at large values of 
n and z and replacing the sum of n contributions by an integral over n. Ultimately, this 
procedure reduces to the following transformations 

1 



J 2 n+1 (z) + Jl.^z) - 2J 2 n (z) -> 



7T 2 e 4 t 2 



a 



dn = — £ dr , 

X 



- > T < OO, 



(A10) 
(All) 
(A12) 
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FIG. 14: The normalized total probability of one-photon emission as a function of £ 2 for different 
values of n max from 25 to 140 as indicated in the legend. Kinematics as in Fig. [T3l 



where <&(?/) is the Airy function, and 



t = (n/2 X ) 2/3 , y = t(l + r 2 ) , X = Z 



kp 



(A13) 




FIG. 15: The normalized total probability of one-photon emission W(£, x) a s a function of x f° r 
different values of £ 2 as indicated in the legend. The crosses and stars correspond to the asymptotic 
values of W(oo,x) for % S> 1 and x ^ 1> respectively. 
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The corresponding probability is expressed in the following form: 



-i/2 



(A14) 



Contrary to Ref. 12| and related papers we do not put £ — > 00 in the above integral, which 
allows to calculate the probabilities at large but finite values of £ 2 in a wide range of x- m 
Fig. [15] we show W(£, x) a s a function of x f° r different values of £ 2 . The crosses and stars 



correspond to the asymptotic values of W(oo, x) at x 3> 1 and x< 1, respectively [12J: 



W(oo,x) = < V , V is 2 y Al5 

^r ( |)(3x 2 / 3 )(i-| r(|)(3 ^ 

In Fig. [16] we show the total probability W(£, x) calculated for a wide region of possible 
experimental conditions ranging from a photon energy of ujl = 0.03 eV up ul = 0.5 keV as a 
function of £ 2 , for £ 2 varying from 0.1 up to 10 4 . This interval covers possible experimental 
conditions illustrated in Fig. CD For £ 2 > 10, W((,, x) is evaluated using the asymptotic 
expression of Eq. (1A14|) . For £ 2 < 10, the probability is calculated as a sum of partial 



harmonics, in our case up to n max = 140. At the matching point £ 2 = 10, the difference 
between the two expressions is less than 5%. 

2. Perturbative treatment 

To prove the equivalence of the probabilities defined above to usual partial cross sections 
we recollect here also the perturbative treatment of the scattering of a photon off an electron 
(Compton scattering, see middle and right panels in Fig. ITT]) . The cross section reads in 
standard notation 

d ^ + '=m^Mf)^ dt - (A16) 

where s = (p + k) 2 is the square of total energy and t = (k — k') 2 denotes the square of 
the momentum transfer. Tfi is the invariant amplitude which consists of direct and crossed 
terms, schematically depicted in Fig. [TT1 (b) and (c), respectively, 

Tfi = e 2 6;(7%(7) [u{p')M^ u{p)\ , (A17) 
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FIG. 16: The normalized total probability of one-photon emission as a function of £ 2 for all 
harmonics summed (solid curves) compared to the summed contributions of the first 25 and 140 
hamonics (dashed curves) with the asymptotic probability given by Eq. (|A14p . The electron energy 
is chosen as E e = 40 MeV, and the photon energies ojl are 0.03 eV (left panel), 1.55 eV (middle 
panel) 0.5 keV (right panel). 

where €(7) and 6(7') denote the polarization vectors of incoming (7)and outgoing (7') pho- 
tons, respectively, u(p) and u(p') are Dirac spinors of incoming and outgoing electrons, 
normalized as uu = 2M e , and M^ u is the transition operator 



' 2p ■ k 11 2p' ■ k ' ' 



(A18) 



In Eq. (IA16I) averaging over the initial and summation over the final spin states are provided. 
The corresponding calculation is described in text books (for example in Q). The result is 
the Klein-Nishina formula 



da 
~dt 



vrr 2 M e 2 



F(p,p',k), 



(s - M 2 ) 2 

where tq = a/M e is the "classical" electron radius and 

Ml 



F(p,p',k) 
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2k -p 2k ■ p' 



Ml Ml 
+ 



2k -p 2k ■ p' 



For the further analysis it is convenient to use the invariant variables u = ^ and dt 



k ■ p k ■ p' 
k ■ p' k ■ p 

k-p 



(A19) 



. (A20) 



2k ■ p du/{l + u) 2 and to employ the probability dW of the 7 + e — > 7' + e' reaction, instead 
of the cross section da, 



dW 



2sp 1 
s + M 2 
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da, 



(A21) 



where p 7 is the photon density (or inverse volume per one photon), defined similarly to the 
dependance of £ 2 on / mentioned in the Introduction, 

s — M 2 M 2 £ 2 

p = _ — 1^^£^. (A22) 

Inspection of Eqs. (1A19|) . (IA21 1) and (1A22[) shows that Compton scattering is independent 
of £ 2 . Variation of £ 2 (or field intensity 7) changes the density of photons p 7 , and therefore 
only determines the overall normalization of dW. 



Appendix B: Perturbative treatment of neutrino pair emission 

Consider first neutrino pair production in the reaction 7 + e — > e' + z^Pj. The differential 
cross section reads 

da® 2 

dM 2 Q du ~ (87r) 3 (s-M e 2 )(l + n) 2 

where Q u is the solid angle of outgoing neutrino in the vv rest frame. The invariant am- 
plitude T"W in lowest order (tree level) of the Glashow-Salam- Weinberg model is described 
by the Feynman diagrams depicted in Fig. 4. Accordingly, the electron neutrino pairs are 
produced via exchange of both the charged W ± and neutral Z° vector bosons, while the 
muon and tau neutrino pairs are produced only through the neutral boson exchange when 
assuming individual lepton number conservation. In the local limit, where all the momenta 
involved in the process are much smaller than the masses of the intermediate vector bosons, 
by making use of a Fierz transformation for the charged boson exchange, one can obtain 
"universal" effective interactions depicted in Fig. [3] as direct (a) and exchange (b) terms. 
This interaction is described by the effective Lagrangian in Eq. fl5]). Then, the invariant 
amplitude is expressed as 

T« = ^ 4%( 7 ) • [u{p')M^ «(p)], (B2) 

where L^a is defined by Eq. (j3J), and the transition operator M^ tMU is defined by Eq. (IA18[) 
with the substitution 

r ^ r{C ®-Cf l5 ). (B3) 



d\l u 

4:71 



(Bl) 
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Summation over spins of the neutrinos (with the assumption m v = 0) and integration over 
the solid angle in Eq. (IB 1 [) leads to 

J ^Tr[L Q Lt] = l{Q a Qp- ga f,Cf), (B4) 
and one can rewrite Eq. (IBlj) as 

dail) ~ ^ ; FU( P ,p',k,Q*) , (B5) 



dM 2 Q du 5127r 2 (s-M e 2 )(l + n) 2 ' 
F®M, k, Q 2 ) = \{g^Q 2 - Q a Q p )g» v ' Tr[M» MfJ]. (B6) 



[1] D. M. Volkov, Z. Phys. 94, 250 (1935) 

[2] V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Quantum Electrodynamics. 2nd ed., 

(Course of theoretical physics; vol. 4), Pergamon Press Ltd. (1982). 
[3] T. Heinzl and A. Ilderton, Opt. Commun. 282, 1879 (2009). 
[4] N. D. Sengupta, Bull. Calcutta Math. Soc. 44, 175, (1952). 
[5] I. I. Goldman, Phys. Lett. 8, 103 (1964) 

[6] L. S. Brown and T. W. B. Klibbe, Phys. Rev. A 133, 705 (1964). 
[7] H. R. Rciss, J. Math. Phys. 3, 59 (1962). 

[8] A. I. Nikishov and V. I. Rhus, Sov. Phys. JETP 19, 529 (1964) [Zh. Eksp. Teor. Fiz. 46, 776 
(1964)]. 

[9] A. I. Nikishov and V. I. Rhus, Sov. Phys. JETP 19, 1191 (1964) [Zh. Eksp. Teor. Fiz. 46, 
1768 (1964)]; 

V. I. Ritus, Sov. Phys. JETP 29, 532 (1969) [Zh. Eksp. Teor. Fiz. 56, 986 (1969)]. 
[10] N. B. Narozhnyi, A. I. Nikishov, and V. I. Ritus, Sov. Phys. JETP 20, 622 (1965) [Zh. Eksp. 

Teor. Fiz. 47, 930 (1964)]. 
[11] A. I. Nikishov and V. I. Ritus, Sov. Phys. JETP 25, 1135 (1967) [Zh. Eksp. Teor. Fiz. 52, 

1707 (1967)]. 

[12] V. I. Ritus, Tr. Fiz. Inst. Akad. Nauk SSSR, 111, 152 (1979) [Proc. P. N. Lebedev Phys. Inst. 

Acad. Sci. USSR 111 (1979)]. 
[13] V. A. Lyulka, Zh. Eksp. Teor. Fiz. 69, 800 (1975). 
[14] N. P. Merenkov, Yad. Fiz. 42, 1484 (1985). 

23 



[15] V. V. Skobelev, Yad. Fiz. 46, 1738 (1987). 

[16] T. Heinzl, A. Ilderton, and M. Marklund, Phys. Rev. D 81, 051902 (2010). 
[17] W. P. Leemans et al, Nature Physics 2, 696 (2006); 

K. Nakamura et al, Phys. Plasmas 14, 056708 (2007). 
[18] A. Henig et al, Phys. Rev. Lett. 102, 095002 (2009). 

[19] J. Faure, C. Rechatin, A. Norlin, A. Lifschitz, Y. Glinec, and V. Malka, Nature 444, 737 
(2006). 

[20] S. P. D. Mangles et al, Nature 431, 535 (2004). 

[21] https : //slacportal . slac . Stanford. edu/ sites/lcls_public/ 

Instruments/SXR/Pages/def ault . aspx 
[22] \protect\vrule widthOpt\protect\href {http : //www. rijnh.nl/f elix/}{http : //www.ri jnh.nl/felix 

[23] G. A. Mourou, T. Tajima, and S. V. Bulanov, Rev. Mod. Phys. 78, 309 (2006); 

Y. I. Salamin, S. X. Hu, K. Z. Hatsagortsyan, and C. H. Keitel, Phys. Rep. 427, 41 (2006); 

M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78, 309 (2006); 

G. Gregori et al, High Energy Density Phys. 6, 166 (2010). 
[24] L. B. Okun, Leptons and Quarks, Amsterdam, Netherlands: North-holland (1982). 
[25] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[26] \protect\vrule widthOpt\protect\href {http : //www . f zd . de/db/ Cms?pNid=132}{http : //www .fzd.de/' 

[27] D. Seipt and B. Kampfer. larXiv:1010.330~T1 [hep-ph]. 

[28] T. Heinzl, D. Seipt, and B. Kampfer, Phys. Rev. A 81, 022125 (2010); 

C. Harvey, T. Heinzl, and A. Ilderton, Phys. Rev. A 79, 063407 (2009). 
[29] N. B. Narozhny, S. S. Bulanov, V. D. Mur and V. S. Popov, Phys. Lett. A 330, 1 (2004). 



24 



